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Unfreezing the radial degree of freedom, we study the three-dimensional Abelian Higgs model with 
compact gauge field and fundamentally charged matter. For small quartic Higgs self coupling and 
finite gauge coupling the model possesses a first order transition from the confined/symmetric phase 
to the deconfined/Higgs phase separated at some hopping parameter k c . Latent heat and surface 
tension are obtained in the first order regime. At larger quartic coupling the first order transition 
ceases to exist, and the behavior becomes similar to that known from the London limit. These obser- 
vations are complemented by a study of the photon propagator in Landau gauge in the two different 
regimes. The problems afflicting the gauge fixing procedure are carefully investigated. We propose 
an improved gauge fixing algorithm which uses a finite subgroup in a preselection/preconditioning 
stage. The computational gain in the expensive confinement region is a speed-up factor around 
10. The propagator in momentum space has a non-zero anomalous dimension in the confined phase 
whereas it vanishes in the Higgs phase. As far as the gauge boson propagator is concerned, we find 
that the radially active Higgs field provides qualitatively no new effect compared to the radially 
frozen Higgs field studied before. 

PACS numbers: 11.15.Ha, ll.10.Wx 



I. INTRODUCTION 



The lattice Abelian Higgs model with compact gauge fields (cAHM) has a wide variety of properties which makes it 
interesting from both the high energy physics [1113 anci condensed matter physics 3] point of view. The compactness 
of the gauge field gives rise to the presence of monopoles which are instanton-like excitations (topological defects) in 
three space-time dimensions. Being in the plasma state, monopoles and antimonopoles guarantee linear confinement 
of electrically charged test particles |j| . The topological defects are forming an oppositely charged double sheet along 
the minimal surface spanned by a Wilson loop. Due to screening, the free energy of this double sheet is proportional 
to the area of the minimal surface. As a result, electrically charged particles experience linear confinement. 

In the case of the cAHM, the plasma state of the monopoles is realized, similar to pure cQED, at small values of 
the hopping parameter (coupling between matter and gauge fields) corresponding to the confined phase of cAHM. 
With respect to the Higgs degrees of freedom that phase is called symmetric. As the hopping parameter increases, the 
system enters the Higgs region, where e.g. the linear potential between test charges is suppressed. Correspondingly, 
in this phase monopoles and antimonopoles arc bound into magnetically neutral dipoles, which provide at best only 
short-ranged interactions between the charged test particles. In other words, the presence of dynamical matter fields 
with charges in the fundamental representation of the gauge group leads to the effect of string breaking which results in 
a flattening of the potential at large charge separations. With respect to the gauge degrees of freedom the Higgs region 
is the deconfined region. The formation of magnetic dipoles in relation to the phenomenon of string breaking were 
studied both analytically in Ref. and numerically in Ref. Q . Note that the origin of monopole binding in the zero 
temperature case of CAHM3 is physically different from the case of monopole binding at the finite temperature phase 
transition occurring in compact QED2+1 without a dynamical Higgs field in the fundamental representation la, El H ■ 

The effect of finite-temperature deconfinement on the photon propagator in (2 + l)-dimensional cQED JgL llCt fll| 
and in zero-temperature cAHM in the London limit are rather similar The propagator is described by a Debye 
mass and by an anomalous dimension which both vanish at the deconfinement transition regardless on its nature, 
which can be caused by finite temperature or by the presence of the matter fields. In both cases the transition is 
attributed to pairing of magnetic monopoles. The mass parameter of the propagator behaves differently in these cases 
since in the deconfined phase of the cAHM the gauge boson acquires a mass due to the Higgs mechanism whereas in 
the case of cQED this mechanism is of course absent. 

The qualitative similarity of the form of the gauge boson propagators in cQED and in the cAHM, so far known in 
the London limit, immediately raises the question about the role of the Higgs field in the emergence of the anomalous 
dimension. Indeed, in the London limit of the cAHM the mass of the Higgs is infinite and the only active ingredient of 
the Higgs field is its compact phase. Away from the London limit the Higgs mass stays finite and the radial component 
of the Higgs field gains influence. In this paper we concentrate on the role of the radially active Higgs field on the 
properties of the gauge boson propagator in conjunction with the phase structure of the model. 
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The phase diagram of the cAHM in three dimension has been extensively studied both perturbatively (analytically) 
and numerically in the literature 0, 0, 0] . Already in the early studies Q the changing nature of the phase 
transition due to radially active Higgs field has been addressed. The authors of Refs. have been concentrated 

on the continuum limit (studying large f3) measuring among others the gauge boson mass from plaquette-plaqucttc 
correlators. As in our previous work we are studying here the propagator in momentum space which allows to establish 
easily whether a non-vanishing anomalous dimension exists. 

The structure of the paper is as follows. Section [H] is devoted to the description of the model and to a sketch of its 
phase structure at two selected values of the coupling describing the quartic self-interaction of the Higgs field. The 
gauge fixing procedure and its inherent problems are reviewed in Section lHfl where also a preconditioning method using 
a finite subgroup of the gauge group is proposed. Near the transition, the changing form of the photon propagator is 
analysed in Section llVl Our conclusions are given in Section Ivl 

II. THE MODEL AND ITS PHASE STRUCTURE 
A. Some properties of the model and its limiting cases 

We consider the three-dimensional Abelian Higgs model with compact gauge fields 9 Xi/i G [— tt,tt) living on links 
I = x,n and with fundamentally charged (Q = 1) Higgs fields <i> x G C on sites x (a vector with integer Cartesian 
coordinates). 

The Higgs field is written in the form 

<!>* = (be*** (1) 
where p x is its radial part and ip x its phase. The model is defined by the action 

S = -P J2 cos 9p ~ K J2 PxPx+fi <*»(-<Px + 0x,m + <At+a) + ^1(pI + KpI ~ !) 2 ) . ( 2 ) 

P X,fJ, X 

where 6p is the plaquette angle representing the curl of the link field 9 x ,n and p- is the unit vector in \i direction. f3 is 
proportional to the inverse gauge coupling squared, j3 = l/(a g 2 ), k is the hopping parameter and A the quartic Higgs 
self coupling. The so called London limit of that model has frozen radial Higgs length corresponding to A — > oo. 

In that limit, at zero value of the hopping parameter k, the model (J2J reduces to the pure compact Abelian gauge 
theory which is known to be confining at any coupling (3 due to the presence of the monopole plasma [4| . We call the 
Iow-k region of the phase diagram the "confined region". At large values of k (also called the "Higgs region") the 
monopoles should disappear because the requirement k — > oo leads to the constraint —tp x + x f + + tp x +(i = 27r/ Xj/i , 
with, in general, l Xttl G ZZ. However, in the unitary gauge, cp x = 0, the above constraint along with the compactness 
condition for the gauge field, 9 X .^ G [— 7T, 7r), gives 9 x ,n = 2ttI x . Pi = 0, which, in turn, indicates the triviality of the 
model at large values of K in the London limit. 

At large (3 the model with A — > oo reduces to the three dimensional XY model which is known to have a second 
order phase transition at k^ y rj 0.453 0] between a symmetric and a Higgs phase. Indeed, in the limit /J^oowe get 
the constraint for the plaquette variable 1 , 9p = [d9] p = 2imp, with np € ffl. The constraint implies dn = (due to 
the nilpotency d 2 =0), which, in turn, gives n = dm, where m £ ^ is a link variable. Thus, the constraint reduces to 
the equation d(6 — 2i:m) = which has a general solution of the form (in usual notations) 9 X4l — — Xx+p, + Xx + 27ri Xl/t , 
where \ x G [— n, it) is a compact scalar field. The integer-valued vector field, l x _^ G Zj, is chosen in such a way that 
9x,n S [— 7r,7r). Fixing the unitary gauge we obtain that in the London limit and at large (3 the action @ reduces to 
the action of the AF-model, S XY = -k^ cos(— \ x + Xx+fi), where the scalar angle x pl a Y s the role of the spin 
field in that model. 

From our previous studies in the London limit at (3 = 2.0 0, Il6j| we know that the transition from the confined 
to the Higgs region is signalled by a rapid drop of the monopole density. In addition, a string breaking phenomenon 
affecting fundamental test charges (being in the same representation as the fundamentally charged matter fields) has 
been observed which are bound by a linear (string) potential in the confined region. While the drop of the monopole 
density signals onset of deconfinement, the anomalous dimension of the photon propagator turning to zero at the 



1 Here and below we often use the differential form notation on the lattice (a, b) = 5Z; a i''i> ll^ll 2 = (fi>@)> the operators d and 5, 
respectively, the lattice curl and divergence. The Laplacian is denoted as A = <5d + d<5. 
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same critical k c shows the transition from the symmetric to the Higgs region, whereas the observed non-zero photon 
mass at larger k is simply due to the onset of normal mass generation by the Higgs effect. 

However, we could not confirm that this crossover is accompanied by an ordinary phase transition. It cannot be 
excluded that the transition is second order with a very small negative critical exponent or possibly of Kosterlitz- 
Thouless type. 

Including a fluctuating radial degree of freedom p x (at small A) changes the situation strongly. For fixed A a first 
order phase transition is known to exist [l3| from the symmetric to the Higgs phase in the /3—k plane accompanied by 
the monopole density approaching zero, similar to our previous studies of the 3D SU(2)-Higgs model |l7j,ll8|. With 
increasing quartic Higgs self coupling the first order transition becomes weaker and ends at a certain critical A c 1191 . 
Above that critical self coupling the phase structure of the compact Abelian Higgs model has been investigated in [IJ 
for relatively large (3 with practically vanishing monopole density. In a very recent study |20j | the location of A c has 
been established for a smaller (inverse) gauge coupling (3 — 1.1 as 

0.030 < A c < 0.032 . (3) 



B. Two regimes of Higgs self coupling 



In the present work we consider the gauge coupling (3 = 2.0 corresponding to a smaller lattice spacing a than in 
Ref. [2(| ■ For this gauge coupling the endpoint of the transition is found to be [2l| 

A c « 0.009(1). (4) 

More details of this phase transition study will be presented elsewhere. We restricted the Higgs self coupling to two 
values A = 0.005 (representing the first order transition regime) and A = 0.02 (representing the continuous transition 
region) . 

In our simulations one complete Monte Carlo update consists of a Higgs and a gauge field update in alternating 
order. For the gauge fields a usual Metropolis algorithm has been applied. The algorithm was adjusted (during the 
thermalization) to ensure an acceptance rate of about 50 per cent. Representing the Higgs field as real two-component 
vector a heat-bath algorithm similar to the method described in j2^ has been used. Additional over relaxation steps 
have been applied to the Higgs field. Since the autocorrelation changes marginally with the number of reflection steps 
not more than two reflections per heatbath update have been used. 

In order to locate the phase transition and to quantify its strength we consider for the present purposes only the 
gauge invariant average Higgs modulus squared 

X 

(L is the linear extent of the lattice) and evaluate its vacuum expectation value (p 2 )- For the Higgs modulus squared 
one can consider the susceptibility 

C P >=L* {({p 2 ) 2 )-(p 2 ) 2 } (6) 

and the Binder fourth cumulant as function of k (at fixed (3 and A) near criticality for different lattice volumes. Widely 
used multihistogram reweighting methods have also been applied in our work in order to determine the pseudocritical 
hopping parameter k c . 

As one example we show here in Fig. ^ for the two A values the histograms P{p 2 ) (normalized to unit area) from all 
data at different K values near k c and lattice sizes L 3 reweighted to the respective pseudocritical hopping parameter 
k c . In case (a) k c for the lower A has been fixed according to the equal height condition applied to the two-state 
histogram. In case (b) for the larger A value k c has been determined according to the condition that the susceptibility 
C p takes its maximum. One observes that in Fig. ^ (a) the two peaks become more pronounced with increasing 
volume whereas the distance between the maxima only marginally decreases. The tunnelling frequency between the 
two coexisting phases in a finite system (obtained from the autocorrelation function of p 2 ) rapidly decreases with 
increasing volume. At a lattice volume 40 3 (the one used for the propagator measurements to be discussed later) no 
tunnelling has been observed. This has allowed us to cleanly study and contrast the respective propagators of the 
competing phases. 

The first order transition is characterized by a finite jump A(p 2 ) of the average Higgs modulus squared, a two-state 
signal at pseudocriticality, surviving in the infinite volume limit. This jump is related to the latent heat of the first 
order transition 

^ = lnPA(p 2 ) (7) 
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FIG. 1: p histograms at pseudocriticality: (a) for A = 0.005 with k c 
A = 0.02 according to the maximal susceptibility criterion. 



fixed according to the equal height criterion and (b) for 



with g| = X/(a0), a is the lattice spacing. In the symmetric phase (p 2 ) is nearly constant as a function of k. In 
contrast to this, on the Higgs side of the transition (p 2 ) increases approximately linear with k. 

Fig. n (a) clearly demonstrates that the transition is first order at A = 0.005. We also have observed that the 
maximum of the susceptibility (not shown here) increases linearly with the volume. Following Ref. fl3T ] we can 
estimate the dimensionless surface tension from the shape of the equal height histogram: 
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(8) 



Here 2L 2 is the dimensionless area of two minimal surfaces splitting the periodic lattice into two regions filled with 
the two pure phases being in coexistence; P ma , x {p 2 ) and P m i n (p 2 ) are the maximal/minimal value of the histogram, 
respectively. An extrapolation to the infinite volume limit (taking into account the leading 1/V corrections) leads to 
the following estimates of the latent heat and of the dimensionless surface tension 



L 



heat 



si 



0.320(2), cr = 0.0121(3). 



(9) 



Fig. ^ (b) shows immediately that at A = 0.02 a first order transition does not exist. As indicated in the figure 
the histogram merges into definitely one maximum at larger lattice volume. Therefore there is no non- vanishing 
latent heat. The width of the histogram reduces with increasing volume. Within errors, the susceptibility for the 
largest studied volumes (again not shown) does not change as function of the volume. Fitting the susceptibility at the 
maximum for lattices from L = 32 to L = 48 by the function C p 2 = AL 1 gives 7 = —0.014(40) with % 2 /d.o.f. = 0.9. 
Including a subleading correction we obtain 7 = —0.012(27) for lattices from L = 16 to L = 48. This excludes a 
second order transition with a critical exponent differing significantly from zero. The transition is similar as the one 
found in the London limit 0- 



C. The monopole density near criticality and autocorrelation estimates 

In addition to p 2 we also consider the monopole degrees of freedom constructed out of the gauge fields. Using the 
differential form notation we consider the monopole charges as the 0-form *j G Z, defined on the sites dual to the 
lattice cubes c [2^| where it is evaluated as follows: 

j = i E (- 1 ) p [M m o d27r - (10) 

Pedc 

The factor (— 1) P takes the plaquette orientations relative to the (outer) boundary of the cube into account, and the 
notation [• ■ ■ ] means taking the integer part modulo 2-7r. The 2-form 



p[j] - m mod2 



(ii) 



5 



corresponds to the Dirac strings living on the links of the dual lattice, which are dual to the plaquettes where the r.h.s. 
is evaluated. The Dirac strings are either forming closed loops or they are connecting monopoles with anti-monopoles, 
<5*p[j] = *j. Whereas *j is gauge invariant, the 2-formp[j] is not. 

The simplest quantity describing the behavior of the monopoles is the average monopole density 



1 



Pmon — T o 



L 3 ' 



(12) 



where j c is the integer valued monopole charge inside the cube c defined in Eq. (|10|l . 

The measured value of the monopole density at smaller and larger A are shown in Fig. [21 for two lattice volumes 
20 3 and 40 3 as a function of k. The vertical lines are drawn to mark the pseudocritical k c . Clearly the monopole 
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FIG. 2: Monopole density p mon near k c for 20 3 and 40 3 lattice at A = 0.005 (a) and A = 0.02 (b). 

density jumps from a nearly constant, finite value to a small value crossing the transition from the confined to the 
Higgs phase. At larger lattices a metastability region is observed (in the form of a small hysteresis cycle) continuing 
the confined and Higgs phase, respectively, into the critical region. For the lattice 40 3 the two metastable states at n c 
are selected by the starting conditions of the Monte Carlo run. At higher A we observe a rapid but smooth decrease 
vs. k near k c precisely as in the London limit. 

Finite size effects on the monopole density are very small, therefore, also the susceptibility of the monopole density 
is practically unchanged with increasing lattice size. As we have found earlier, the remaining non-vanishing density 
above the transition is due to dipole pairs forming a dilute dipole gas. We will keep to these two values of the quartic 
self coupling in order to study the influence of the different transition regimes on the photon propagator. 

Finally we show in Fig. [2Jtlie integrated autocorrelation time to (for p 2 and p mon ) as function of k. We recall that 
the autocorrelation time of an observable O is defined as 



1 - 

ro = ^ + ^C (r), AT>! (13) 



in terms of the autocorrelation function Co(t). In general, Co{t) is defined as 

(O t O t+T ) (O t ) 2 

Go{t) = m <o t ) 2 (14) 

where the observable O as measured at Monte Carlo "time" t is folded with O measured at t + r. On the confinement 
side the autocorrelation time for the monopole density is much larger than for the Higgs modulus squared. For both 
quantities to reaches 0(100) in the first order regime and up to 0(1000) at critical coupling at larger A. Since for 
the latter case the autocorrelation times of both L = 20 and L = 40 cubic lattices are roughly the same, we have 
not to worry about critical slowing-down. These autocorrelation times have been used to select statistically (almost) 
independent configurations for measurements of the photon propagator in the respective cases of the quartic Higgs 
self couplings. 
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FIG. 3: The integrated autocorrelation times to of p 2 (circles) and p mon (squares), (a) in the first order region, (b) in the 
continuous transition region. 



III. LANDAU GAUGE FIXING WITH OVERRELAXATION REVISITED 
A. Landau gauge and Dirac strings 



In order to calculate the (gauge dependent) photon propagator directly, a special gauge has to be chosen, in our 
case the minimal Landau gauge 2 . This gauge is defined by finding the global maximum of the gauge functional 



F 



1 



' Minks ^ 



• max . 



(15) 



The gauge field and Higgs field angles have to been transformed into this gauge within gauge transformations G 
with elements g x = exp(iu> x ) acting as 



fx = (fx + u x ) 27T 



(16) 
(17) 



until the condition (|15|l is satisfied. The notation (. . . )„ indicates that values are restricted to the interval [— 7r,7r) 
(modulo 2ir). 

Every gauge transformation can be decomposed into a sequence of local gauge transformations on each lattice site 
which would increase the sum l|15|) restricted to neighboring links. The corresponding angle oj x is found easily as 



tanu} x = ■ 



(18) 



The efficiency of the full optimization procedure depends crucially on the landscape of the gauge functional in the 
space of all gauge transformations. We have applied the local gauge transformation in a checkerboard fashion (based 
on a separation into odd and even sublattices). It turns out that overrelaxation allows us to speed up the finding of a 
local maximum with respect to full gauge transformations G. This was realized multiplying the gauge angles uj x by 
an overrelaxation parameter 1 < rj < 2. Note that the extreme case rj = 2 does not change the local gauge functional 
at all because it corresponds to a microcanonical update with respect to the gauge functional F. Fastest convergence 
was obtained for r\ values of about 1.8-1.9 in agreement with earlier findings jlfj. 

This algorithm will in general not lead to the absolute (global) maximum of the gauge functional F l|15|) . Typically 
it will get stuck in one of the local maxima, which are called Gribov copies of each other and of the (unknown) 



2 For more details on the Landau gauge see llOl and references therein. 
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true maximum. As long as this non-uniqueness influences the non-gauge-invariant observable of interest, this is the 
origin of the so-called Gribov problem. Practically, it is partially cured by repeating the same gauge fixing procedure, 
applying it to random gauge copies of the original Monte Carlo configuration. The local maximum with the largest 
value of F is taken as the tentative global maximum. This rests on the assumption that one of these random gauge 
copies might be situated in the basin of attraction of the global maximum. 

The number of gauge equivalent configurations produced to restart the gauge fixing is denoted as Ng, and the 
iterative gauge fixing generically leads to really different maxima. We have then to content with the best out of all 
Nq + 1 local maxima of the gauge functional, and the observable of interest (the photon propagator) will depend on 
N G . 

For the purpose of the following discussion 77 = 1.8 has been chosen. It will become clear that this choice is not 
advantageous. A closer investigation of the gauge fixing procedure reveals a direct relationship between the achieved 
value of the gauge functional and the density of leftover (gauge dependent) Dirac strings 3 



iV links 



(19) 



where the sum is taken over the links * P of the dual lattice. In particular, maximizing the gauge functional leads to a 
minimization of the Dirac string density as discussed in |10| . Note that this dependence is a result of the local gauge 
fixing procedure which blindly acts with respect to the Dirac string content. 

The different behavior in the confined and Higgs (deconfined) phase is demonstrated in Fig. 0] showing the Dirac 
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FIG. 4: (a) The highest gauge functional after 1 trial (Ng = 0) and 171 trials (Ng = 170) on the confined side of the continuous 
transition near k c ; (b) the same with 1 and 61 trials, respectively, on the Higgs (deconfined) side. 



string density (separately measured with respect to different directions, but shown only for the x direction) vs. the 
values of the gauge functional for an ensemble of 195 independent configurations. Both quantities have been finally 
achieved after inspecting the first or the best copy among Nq + 1 local maxima of the gauge functional F. 

From the figures it is clear, that in both regions with a single attempt the best realization of the gauge fixing has 
never been found. In the Higgs region corresponding to deconfinement (k = 0.378), with sufficiently many copies 
Nq created for selection, in most cases the relative maximum of the gauge functional is located in a small interval 
just below 0.955 with a Dirac string density below 0.001. Here one can hope finally to catch the global maximum 
accompanied by a minimal density of Dirac strings with a modest increase of the number of trials. 

In the confining phase the behavior is different. Even with 170 copies under inspection, the reduction of the Dirac 
string density and the increase of the gauge functional is less efficient. Although the average of the gauge functional is 
shifted towards higher values its spread compared to the original spread remains rather large and the highest values 
of 0.938 are scarcely reached. The relative decrease of the density of Dirac strings is poor. We will see soon that such 
a far-from-perfect gauge fixing strongly influences the propagator. 



3 For the definition of the Dirac strings see Eq. 1111 . In four dimensional compact U(l) the attention was drawn to Dirac sheets in |24| . 
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On the other hand, from the absolute values of noirac it can be concluded that Dirac lines wrapping around the 
torus are absent. They cannot be blamed to spoil the photon propagator, an effect which otherwise could be argued to 
be unphysical 4 . Essentially the same observations can be made in the case of the smaller A value (in the neighborhood 
of the first order transition), but here without the possibility to simply separate out the "physical case". Thus, it is 
conceivable that a smarter local gauge fixing procedure could be able to reduce the Dirac string density to a (gauge 
independent) minimum. 

Before any gauge fixing is attempted, the gauge dependent Dirac strings are either open lines (forming Dirac lines 
connecting monopoles and antimonopoles) or closed loops. In the confined phase these strings form huge clusters 
where the individual lines and loops cannot be resolved. It is impossible to uniquely identify throughgoing lines/loops 
at points where more than two Dirac strings meet. After some local maximum of the gauge functional has been 
found, the number of closed loops has decreased considerably. Closed loops have the chance even to be gauged away. 
Open strings which connect gauge-independent monopolc-antimonopole pair positions are only able to minimize their 
length (number of contributing strings) during the gauge fixing procedure. 

In an idealized picture with no leftover loops, a local maximum of F is associated with a realization of Dirac 
lines necessary to connect monopole pairs. The global maximum then corresponds to a minimal realization. In the 
confined phase we have a dense plasma of monopoles and antimonopoles with total charge zero. Hence a lot of 
pairings are possible which corresponds to a huge number of local maxima of the gauge functional. On the Higgs side 
(deconfinement) only a small number of nearby (in lattice distances) monopole-antimonopole pairs (a dilute dipolc 
gas) is left. Consequently, the number of local maxima is much lower than in the confined case. This picture explains 
why it is much more unlikely to find the global maximum of F in the confined phase compared to the deconfincd 
phase. The problem of the gauge fixing in the confined phase belongs to the class of complex systems similar in 
complexity to the zero temperature states in spin glasses. 

B. Choice of the algorithm 

The gauge fixing algorithm with continuous compact fields is very costly. To speed up the algorithm several 
refinements have been done. The implementations were tested for n < k c (confined side) because this is the most 
critical region. Firstly, we tested the effect of restricting the link variables to the discrete subgroup Zjv C U(l), 
i.e. mapping first the gauge field variables to the closest Zaj- realization. Consequently, the gauge transformations 
were restricted to Zat, too. We observed that approximating gauge field angles by I byte arrays (N = 256) gave, on 
the average, a lower gauge functional than using continuous compact fields. Using 2 byte arrays (N = 65536) the 
gauge functionals resulting from gauge fixing using continuous angles and discrete angles, respectively, coincide within 
statistical errors for some tested rj values. 

Secondly, we used a preselection strategy. The Nq + I gauge fixing attempts can be performed at some chosen 
overrelaxation parameter r\ < 2 with a rather weak stopping criterion, i.e. the lower limit for the change of the 
gauge functional was still relatively large. The actual choice of the best r\ for the preselection stage is discussed below. 
Afterwards, among the non-precisely gauge-fixed copies the trial with the highest gauge functional was taken up again, 
and a final gauge fixing with a strong stopping criterion was applied stepping back now to a smaller overrelaxation 
parameter r\ — 1.90 leading to fastest final convergence during the final-gauge fixing stage. The intuition behind 
this method is that during the preselection stage a certain pairing of (anti-)monopoles is already chosen and the fine 
tuning of the link angles is left to the second stage. 

Both ideas have been then combined to use discrete gauge fixing in the preselection stage and to apply continuum 
gauge fixing only to the best copy obtained in the result of preselection. For this purpose, a Zat realization of the 
gauge field was used first to construct a suitable Zat gauge transformation, applying the weak convergence criterion. 
Then, embedding the optimal discrete gauge transformation into U(X) as an initial guess and returning to the true 
U(l) gauge field, the final continuous gauge fixing was done only for the most promising copy. The described combined 
procedure is called optimized gauge fixing in Fig. [3] Sensible stopping limits depend on the gauge functional landscape 
and thus on the lattice size and the coupling parameters. They had to be chosen appropriately. Using that optimized 
gauge fixing a speed-up of roughly a factor 10 has been found. 

As mentioned above using rj = 1.8 might be not a good choice. This is demonstrated in Fig. [5] showing the rj 
dependence of the gauge functional F and of the Dirac string density riDirac for 77 values near (but below) two using 
20 independent measurements. This dependence has been overlooked in our previous studies. Only for the largest 77 



4 See Ref. Illl on the photon propagator in the deconfined phase of CQED3 where it was possible do deal with this unphysical remnant 
of imperfect gauge fixing. 
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FIG. 5: (a) The average of the maximal gauge functional per configuration as function of 77; (b) the same for the minimal 
Dirac string density. Continuous gauge fixing is selfexplaining. For details of the optimized gauge fixing see the text. 



values (very close to 77 = 2) we seem to reach a new level of performance (higher gauge functional, lower Dirac string 
density). In the figure for some selected r\ values we also show measurements comparing the results achieved with 
the optimized algorithm (using discrete stage for preselection and preconditioning) with the standard algorithm using 
only continuous gauge transformations. As mentioned above, within statistics the results coincide. At 2 — 77 = 0.02 we 
show in addition how the change of Nq influences both gauge functional and Dirac string density. In order to reach 
the same level by increasing the number of copies Nq turned out to be much more expensive in computing time than 
increasing 77. Choosing 77 close to 2 or, alternatively, applying several 77 = 2 (the exactly microcanonical case) and 
77 < 2 overrelaxation steps in alternating order, has led to a very similar convergence behavior. In the first case the 77 
value was tuned and in the latter case we have adjusted the number of microcanonical steps interspersed between the 
overrelaxation steps in order to obtain the largest average of the final gauge functional. The iteration number 7V iter 
(iteration of local gauge transformations) needed in order to reach a certain gauge functional has been found to be 
about the same in both cases. 

Unfortunately, iVjter increases very fast if one tries to get at higher values of the gauge functional using one of 
these methods. Comparing 77 = 1.997 with 77 = 1.80 that number rises by a factor of 15 to 20. The highest value 
of 77 which has led to an affordable Niter in order to complete at least 50 independent measurements for the photon 
propagator was actually 77 = 1.997 . Using this value one reaches the edge of a new level of performance both in the 
gauge functional and the Dirac string density indicating that the effect of increasing 77 further should become small. 
Note that for 77 close to 2 the algorithm behaves similar to a simulated annealing technique. Which of the techniques 
is computationally favorable, remains to be studied. It has been observed that increasing 77 is much more efficient 
than increasing the number of trials Nq + 1 using as measure the computational workload, i.e. the total number of 
updates of the local gauge transformations (Nq + 1) x Mtcr- 

After these detailed studies we have chosen the following parameters as a standard to perform the Landau gauge 
fixing in the measurements of the gauge boson propagator: As weak stopping criterion the change of the gauge 
functional between two iterations had to become less than 1/(3 x 40 3 ), as strong stopping criterion we used 10 _8 /(3 x 
40 3 ). The number of Gribov copies at 77 = 1.997 was chosen to be 51. 



IV. RESULTS FOR THE PROPAGATOR NEAR CRITIC ALITY 



The photon propagator is the gauge-fixed ensemble average of the following bilinear in A, 

D flv (f) = (A^A_^), 



(20) 
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where A is the Fourier transformed lattice gauge potential (in lattice momentum space) related to the gauge links 
9 X u (fixed to minimal Landau gauge in coordinate space) via 5 

h» = ( j3) V2 E ex P( 2 ^ E ^ } ) Sm ^' *m = 0,±1,...,±^, (21) 

^ ' X v—\ 

The vectors of lattice momenta p on the left hand side of H20|) are related to the integer valued Fourier momenta 
vector k = (fci, &2, fe) as follows: 

2 irk 

P P .(k fJ ,) = -sm—^-, p=(p 1 ,p 2 ,p 3 ). (22) 
a Li^ 

Thus the lattice equivalent of the same continuum p 2 can be realized by different vectors k which eventually could 
reveal a breaking of rotational invariance. 

Assuming reality and rotational invariance the most general tensor structure of the continuum propagator is 

D,AP) = D( P 2 ) + P ^ F(P P (23) 

with the three-dimensional transverse projection operator 

P»u(p)=5 fiU - P ^. (24) 

The two scalar functions (form factors) D(p 2 ) and F{p 2 ) can be extracted on the lattice from D^ y {p) by projection. 
Since D^ v (p) is only approximately rotationally invariant, the form factors may be scattered rather than forming a 
smooth function of p 2 . When the Landau gauge is exactly fulfilled F(p 2 ) = 0. On the lattice, this is actually the case 
as soon as one of the local maxima of the gauge functional (|15Jl is reached, with an accuracy which directly reflects 
the precision at stopping of the gauge fixing iterations. There is no possibility to monitor through F(p 2 ) to what 
extent the global maxima have been successfully found. 

In order to measure the propagator, for lattice size 40 3 , we have considered 50 (almost) independent configurations 
(separated by 720 Monte Carlo sweeps) on the confined/symmetric side and 100 configurations (separated by 360 
Monte Carlo sweeps) on the deconfined/Higgs side near k c for the two representative cases of A (A = 0.005 and 
0.02). There were several reasons for choosing this lattice size. A large lattice reduces the influences of the boundary 
and leads to smaller variances in certain observables, e.g. the monopole density. Furthermore, in the first order 
regime tunnelling was prohibited (mixed phases during measurement) and the propagator properties of both phases 
in a metastability region could be independently studied. Finally, the relatively large size allows to measure the 
propagator for more realizations of very small momenta where the propagator in the confined phase is most sensitive. 

In FigElwe present as example the form factor D(p 2 ) multiplied by p 2 measured at the same k = k c in the metasta- 
bility region. The system was found either in the confined/symmetric and the deconfined/Higgs phase depending on 
the initial conditions of the Monte Carlo run. 

To describe the scalar function D(p 2 ) quantitatively we applied the following fit ansatz 

g ui 2a 

D{P2) ^ /3 [p2(l +Q ) +TO 2(l +a) ] +C ' ( 25 ) 

where Z : a, m and C are the fitting parameters. This fit has been successfully used to describe the propagators 
of the finite and zero-temperature compact U(l) gauge mo del 01, ITof in 2 + 1 or 3 dimensions, respectively and the 
propagator in the London limit of the Abelian Higgs model |l6j|. The form is similar to some of Refs. [2{| where the 
propagator in gluodynamics has been studied. The meaning of the fitting parameters in Eq. (|25|) is as follows: Z is 
the renormalization of the photon wavefunction, a is the anomalous dimension, m is a mass parameter. As shown 
in Ref. ^(Jj m CQED3 this mass parameter coincides with the Polyakov prediction for the Debye mass, generated 
by the monopole- antimonopole plasma. The parameter C corresponds to a (5-like interaction in the coordinate space 
and, consequently, is irrelevant for long-range physics. We should remark that before fitting the data, the propagator 
values have been averaged over all those lattice momenta realizing the same p 2 . 



The propagator in the continuum formulation is given by (/3/ 'a)D ^(p). 
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(a) (b) 

FIG. 6: Momentum dependence of the photon propagator form factor D(p 2 ) multiplied by p 2 for A = 0.005 at re c — 0.3579 
together with the fit l!2.~>> in the confined/symmetric phase (a) and in the deconfined/Higgs phase (b). 

The momentum dependence of the gauge boson propagator is well described by the fitting function both in the 
confined/symmetric and in the deconfined/Higgs phase. The enhancement at intermediate momenta disappears 
(abruptly or continuously) with increasing re. In contrast, the suppression at very small |pj remains at all re values. 
We can conclude that the photon propagator is always less singular in the infrared than a free one. According to our 
fitting function (|25(l the propagator is finite at p 2 = 0. 

In Figs. we demonstrate how the momentum dependence of the propagator changes when crossing the first order 




(a) (b) 

FIG. 7: Momentum dependence of the fitted photon propagator form factor D(p 2 ) multiplied by p 2 for different re's near 
criticality for the first order (a) and continuous transition (b) regime. 

or continuous transition. In the confinement region sufficiently away from the critical region the propagator form 
factor is practically the same irrespectively of the A-value. Approaching the critical n c from below, the form factor 
p 2 D(p 2 ) remains enhanced at intermediate momenta in the first order regime of small A, changing then abruptly to 
a free massive propagator in the Higgs phase (deconfincment). Contrary to that, p 2 D(p 2 ) varies continuously with 
increasing n crossing the (continuous) transition at the higher A value. 

In order to demonstrate the influence of a less efficient gauge fixing procedure on the shape of the emerging photon 
propagator we show in Fig. [8] the rj dependence of the fit parameters Z and a on the confined side in the case of 
larger quartic Higgs self coupling obtained in 20 measurements. The dependence of the fit parameters resembles that 
of the gauge functional and the Dirac string density presented in Fig. [3] with a transition from a less "perfect" level 
of performance to a better one around 2 — ?y = 10 _1 ...10 -2 . This strongly suggests that these fit parameters are 
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FIG. 8: The fitted wave function normalization Z (a) and anomalous dimension a (b) at A = 0.02 and k = 0.376 as functions 
of the overrelaxation parameter rj. 



closely related to the Dirac string density. We have checked that the minor discrepancy between the fit parameters 
at r) — 1.85 using the optimized and continuous gauge fixing goes away with increasing statistics. 

In contrast to this, the other fit parameters to and C are roughly rj independent which points out that they depend 
only on the (gauge independent) monopole content. Still, a minor r\ dependence of Z and a might remain because 
we have used r\ = 1.997 as a "still affordable" choice. Probably a different gauge fixing algorithm has to be used to 
reduce this influence below the one-percent level. Comparing the obtained fitting results with those obtained earlier, 
with overrelaxation parameter rj = 1.8, for pure U(l) and for the Abelian Higgs model in the London limit 0] 

we have to admit that the presented numerical values for Z and a in the confined phase are likely to reach values 
smaller by roughly 20 to 25 per cent once a better tuned rj is used. However, all qualitative results remain unchanged. 

The resulting fit parameters near the first order transition are presented in Fig. [5] as functions of k. All parameters 
show a discontinuity at the phase transition. The metastability can be seen clearly reaching the transition from 
below (confined/symmetric phase) or from above (deconfined/Higgs phase). The anomalous dimension jumps from 
a non-zero positive value at low n to zero in the Higgs phase. The increasing mass of the photon propagator in the 
Higgs phase arises from the Higgs mechanism, whereas in the case of confinement the non-zero mass is due to the 
monopole plasma. Note that the mass is only very weakly dependent on k up to n c in agreement with the behavior 
of the monopole density [compare Fig.[2Ja)]. 

The corresponding results for the continuous transition regime are presented in Fig.^| The fit parameters behave 
in complete analogy to the London limit: In the region very close to k c , the mass to shows a minimum caused, as 
one could guess, by the interference between perturbative mass and Debye mass effects. Indeed, as k tends to k c , the 
Debye mass gets smaller since the density of the monopole-antimonopole plasma drops rapidly. One the other hand, 
the perturbative mass term becomes more significant. The interplay of these two tendencies results in the noticed 
minimum at K Rs k c . 

We have fitted the behavior of the mass parameter, to, as function of k on the Higgs side near the transition by 
the function 

to(k) = h(k — k^*) , k > k c . (26) 

The fit ansatz implies a vanishing mass at k c what is actually not found numerically. So we treat that ansatz only 
phenomenologically. The best fit with x 2 /d.o.f. — 1.3 allows to determine = 0.3767(3), which agrees with the value 
from thermodynamic measurements given in Fig.^b) within errors. The estimated "critical exponent", v = 0.19(6), 
differs from the mean field exponent, vmf — 1/2 which had been confirmed in Ref. 0, 0] for larger (3. The fit 
including the four nearest fitted masses is shown in Fig. HOf c) by the solid line. 

From the results for the anomalous dimension in both regimes we conclude that in cAHM with radial degrees of 
freedom the effect of the monopole pairing on the propagator is the same as in cQED: the anomalous dimension gets 
close to zero in the Higgs phase which is dominated by the magnetic dipole gas. This change is continuous in the 
continuous transition case (large A), but discontinuous for the first order phase transition (small A). 

Similarly to the anomalous dimension, the effect of the monopole pairing on the renormalization parameter Z of the 
photon wavefunction in cAHM, shown in Fig. llOf a). is remarkably similar to the cQED case observed in Refs. |9HTo|. 
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FIG. 9: Fit parameters for the propagator at f3 = 2.0 and A = 0.005 (in the first order transition region) as function of k: (a) 
Z; (b) a; (c) m; (d) C. 



The wave function renormalisation Z drops continuously across the continuous transition, but discontinuously at the 
first order transition for k reaching k c from below. 



V. CONCLUSIONS 



We have carefully investigated the problems afflicting the Landau gauge fixing procedure based on overrelax- 
ation. An improved gauge fixing algorithm has been proposed and tested which uses a finite subgroup in a preselec- 
tion/preconditioning stage. The computational gain in the expensive confinement region is a speed-up factor around 
ten. Such a stage could be useful for non-Abelian gauge groups as well. 

We found that the gauge boson propagator of the three-dimensional compact Abelian Higgs model with radially 
active Higgs fields possesses a non-zero positive anomalous dimension in the confinement phase (at small k). The 
momentum dependence of the propagator can be described by four parameters: the mass, the anomalous dimension, 
the photon renormalization function and the strength of the contact term. 

At the first order phase transition (small Higgs self coupling) the parameters of the propagator show a discontinuity 
while in the vicinity of the continuous transition (large Higgs self coupling) the change of the parameters is continuous. 
Apart from the difference caused by the nature of the phase transition, the propagator behaves qualitatively similar 
to the propagator calculated in the London limit. Thus we conclude that the main ingredient which influences the 
gauge boson propagator is played by the compact phase of the Higgs field. 

Note that the anomalous dimension of the gauge field does not become clearly negative neither in the confined 
phase nor in the deconfined phase. As it is found in the case of the cQED [jj, Qij the compactness of the gauge field 
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provides a positive contribution to the anomalous dimension while the role |3|, |26J of the dynamical Higgs fields is to 
make the anomalous dimension negative. Our study suggests that the positive contribution due to the compactness 
of the gauge fields is stronger than the negative contribution from the Higgs fields. 
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